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\  /  ABSTRACT 

'  I 

*<ll 

Experimentally  it  had  been  observed  in  1969  for  cyclic  ethers  with  con¬ 
ventional  substituents  that  their  propensity  to  polymerize  was-  influenced 
by  their  basicities  and  ring  strain.  For  energetic  new  polymers  it  was 
desired  to  predict  the  propensity  to  polymerize  of  cyclic  ethers  substituted 
with  energetic  groups  prior  even  to  synthesis  of  the  monomers  themselves. 

We  carried  out  ab-initio  quantum  chemical  calculations  on  these  energetic 
cyclic  ethers  and  generated  the  electrostatic  molecular  potential  contour 
maps.  These  maps  indicate  vividly  the  basicities  of  the  cyclic  ethers  and 
permit  the  relative  ranking  of  their  tendencies  to  polymerize.  The  maps  also 
allow  optimal  choice  of  good  copolymer  candidates. 
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I.  Introduction 

In  support  of  an  experimental  study ^  concerned  with  the  synthesis 
and  polymerization  of  monomers  of  cyclic  ethers  substituted  with  exotic 
energetic  groups  such  as  nitro,  fluoro,  azido  and  azidomethyl  we  have  carried 
out  quantum  chemical  calculations  designed  to  give  insight  on  the  following 
even  prior  to  the  synthesis  of  the  monomers  themselves. 

1)  Propensity  of  cyclic  ethers  to  undergo  cationic  polymerization 

2)  What  happens  as  a  catalyst  (a  Lewis  acid)  approaches  the  cyclic 
ether  monomer  and  what  are  the  preferable  catalysts. 

3)  What  happens  when  a  complex  of  cyclic  ether  plus  catalyst  approaches 
another  molecule  of  cyclic  ether  monomer. 

4)  In  some  instances  polymers  derived  from  single  monomers  yield  a 
product  that  is  too  refractory.  In  such  an  event  what  would  be 
suitable  and  optimal  reactants  for  co-polymerization  to  avoid 
this  undesirable  property. 
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It  had  been  observed  by  Japanese  investigators  who  studied  experimen¬ 
tally  three-to-six  membered  cyclic  ethers  with  alkyl,  alkoxy  and  aryl  sub¬ 
stituents  that  the  rate  of  polymerization  of  cyclic  ethers  was  proportional 
to 

a)  the  base  strength  of  the  cyclic  ether 

b)  the  ring  strain  of  the  cyclic  ether 

Thus,  to  determine  propensity  for  polymerization  it  was  necessary  to  find 

theoretical  indices  that  would  be  predictive  of  basicity  and  ring  strain. 

3  4  5 

The  studies  of  Scrocco  and  Tomas i  *  ,  Petrongolo  and  other  investi¬ 
gators®  on  the  electrostatic  molecular  potential  contour  maps  around 
molecules  have  demonstrated  successful  application  of  such  contour  maps  to 
a  wide  variety  of  chemical  problems  (reactivity,  structure,  scattering  and 
energetics  of  organic,  inorganic  and  biological  systems). 

Our  previous  ab-initio  quantum  chemical  calculations  on  morphine  (an 
N-methyl  narcotic  agonist)  and  on  nalorphine  (an  N-allyl  narcotic  angonist- 
antagonist)^  had  shown  that  while  the  N  atom  in  both  these  compounds  had 
the  same  charge  (which  we  also  verified  experimentally  by  measuring  their 

ESCA  spectra®)  morphine  and  nalorphine  had  very  different  pKfl's  and  also 
15  9 

different  N  NMR  shifts.  The  electrostatic  molecular  potential  contour 
maps  we  generated  from  the  ab-initio  wave  functions  around  the  N  atoms  in 
morphine  and  nalorphine^®  showed  that  while  the  charge  on  the  N  atom  was 
the  same  and  the  maps  were  similar  in  the  N  lone  pair  direction,  the  maps 
around  the  N  atoms  were  quite  different  in  other  directions.  Obviously, 
charge  alone  on  an  atom  does  not  seem  to  be  sufficiently  sensitive  to 
distinguish  different  degrees  of  basicity,  especially  for  closely  related 
molecules.  However,  the  use  of  electrostatic  molecular  potential  contour 
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maps  did  hold  promise  of  being  a  sensitive  criterion  for  basicity.  Comparisons 
of  the  sizes  of  the  negative  regions  of  the  maps  around  the  0  atom  in  the 
ring  should  be  indicative  of  the  relative  basicities  of  these  exotic  energetic 
cyclic  ether  monomers.  The  3-  and  4-  member  rings  are  the  ones  of  choice 
for  the  experimental  program.  To  this  end  we  have  carried  out  ab-initio 
MODPOT/VRDDO/MERGE  calculations  on  oxetane,  3-fluoro-3-nitrooxetane  and 
3,3-dinitrooxetane  (which  were  of  first  interest  experimentally)  and  gener¬ 
ated  the  electrostatic  molecular  potential  maps  around  the  molecules.  Maps 
were  also  generated  around  3-azidooxetane  (of  next  interest  experimentally) 
and  around  the  complex  of  BF^  catalyst  attacking  3-azidooxetane.  The  results 
and  interpretations  of  these  studies  are  reported  here. 
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II.  Computational  Details 

11  12 

The  computations  were  carried  out  with  our  recent  ab-initio  programs  ’ 

13-15 

that  also  incorporate  several  desirable  options  :  ab-initio  effective 

core  model  potentials  (MODPOT)'®,  which  enable  one  to  treat  the  valence 

electrons  only  explicitly,  yet  accurately;  a  charge-conserving  integral 

prescreening  evaluation^  (which  we  named  VRDDO-variable  retention  of 

diatomic  differential  overlap*^  especially  effective  for  spatially 

extended  molecules;  and  an  efficient  MERGE  technique  to  reuse  integrals 

14e 

from  common  skeletal  fragments  . 

The  ab-initio  effective  core  model  potential  formalism  used  is  the 
one  introduced  by  Bonifacic  and  Huzinaga^.  The  model  potential  consists 
of  a  static  local  potential  that  simulates  the  effect  of  the  Coulomb  and 
the  exchange  of  the  atomic  core  electrons  and  a  set  of  projection  oper¬ 
ators  that  prevent  the  collapse  of  the  valence  shell  orbitals  into  the 
inner  cores.  Only  the  valence  electrons  are  considered  explicitly.  In 
the  molecular  case  they  introduced  a  spherical  static  potential  to  replace 
the  local  potential. 

Our  previous  experience  had  shown  that  ab-initio  MODPOT/VRDDO  results  would 

agree  to  about  three  decimal  places  for  valence  orbital  energies  and 

population  analyses  and  to  better  than  three  decimal  places  for  isomer 

total  energy  difference  and  potential  energy  curves  with  those  quantities 

calculated  by  completely  ab-initio  calculations  using  the  same  valence 
14 

atomic  basis  set  .  The  MERGE  technique  enables  considerable  savings  in 
computer  time  when  one  is  studying  a  series  of  similar  compounds  of 
even  medium  size. 
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The  definition  of  the  electrostatic  potential  of  a  molecule  is 
given  by 

nuclei  +  +  -►  + 

V(r)  =  E  za/lra  "  rl  -  /dr1p(r1)/|r1  -  r| 

a 

where  p(r-j)  is  the  electron  distribution  derived  from  the  molecular  wave 
function  and  are  the  nuclear  charges.  The  first  order  interaction  energy 
between  the  molecular  charge  distribution  and  a  point  charge  distribution 
is 

W  =  Z1q(ri )V(r7 ) 

Such  an  expression  has  previously  been  used  for  comparative  purposes. 
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III.  Results  and  Discussion 

Ab-initio  MODPOT/VRDDO  and  ab-initio  MODPOT / VRDDO/ ME RGE  calculations 
were  carried  out  for  oxetane,  3-fluoro-3-nitrooxetane  and  3,3-dinitroxetane 
(Figure  1),  the  geometries  were  optimized  and  electrostatic  molecular 
potential  maps  were  generated  around  the  molecules  in  three  directions: 
i)  parallel  to  the  mean  molecular  plane  (xz);  ii)  parallel  to  the  (xy) 
plane  and  iii)  parallel  to  the  (yz)  plane.  While  the  potentials  were 
calculated  for  a  1000  points  per  plane  for  each  of  several  dozen  planes 
around  the  molecule  we  have  only  shown  the  maps  of  two  of  the  most  illustra¬ 
tive.  Also  on  these  maps  the  contours  were  drawn  down  to  the  -20  kcal 
contour  (rather  than  include  the  more  negative  contours)  so  as  not  to 
clutter  the  maps.  By  focussing  on  the  -20  kcal  contour  and  watching  the 
size  of  that  contour  it  is  illustrative  of  the  relative  difference  in 
basicities.  The  results  for 

(a)  Oxetane  (Figures  2a  and  b) 

(b)  3-Fluoro-3-nitrooxetane  (Figures  3a  and  b) 

(c)  3,3-Dinitroxetane  (Figures  4a  and  b) 

indicate  quite  clearly  that  oxetane  is  the  most  basic,  3-fluoro-3- 
nitrooxetane  is  next  and  3,3-dinitrooxetane  is  the  least  basic.  The  ring 
strains  are  not  expected  to  differ  nearly  as  drastically  in  these  molecules 
of  similar  ring  size  as  do  the  basicities.  We  are  currently  developing 
theoretical  procedures  for  calculating  or  estimating  relative  ring  strains 
in  these  cyclic  ethers.  Thus,  in  response  to  question  1),  it  would  be 
expected  that  oxetane  would  polymerize  more  easily  than  fluoronitrooxetane 
and  that  more  easily  than  dinitrooxetane.  Our  calculations  were  carried 
out  prior  to  the  experiments  to  polymerize  fluoronitrooxetane  and  dinitro- 
oxetane.  Very  recent  experimental  results  validate  our  predictions  . 
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It  was  already  known  from  the  previous  Japanese  work  that  oxetane  would 
polymerize  with  a  Lewis  acid  catalyst.  In  the  experimental  work  carried 
out  in  the  overall  ONR  project  3-fluoro-3-nitrooxetane  was  synthesized 
and  it  undergoes  cationic  polymerization  with  a  Lewis  acid  catalyst. 
3,3-Dinitrooxetane  was  also  synthesized,  however  its  basicity  is  so 
diminished  that  it  is  extremely  difficult  to  polymerize.  The  polymerization 

of  3 , 3-di ni trooxetane  cannot  be  initiated  in  the  conventional  way  with  a 
Lewis  acid,  no  matter  how  strong  the  Lewis  acid. 

We  are  currently  also  carrying  out  similar  calculations  to  generate 
electrostatic  molecular  potential  contour  maps  around  a  variety  of  other 
exotical ly  energetically  substituted  oxetanes.  These  maps  are  sufficient 
to  indicate  the  propensity  of  these  oxetanes  to  undergo  initiation  for 
cationic  polymerization.  If  the  maps  indicate  a  basicity  of  the  substi¬ 
tuted  oxetane  at  least  as  great  as  that  of  3-fluoro-3-ni trooxetane  the 
monomer  will  undergo  cationic  polymerization.  If  the  map  indicates  the 
basicity  of  the  monomer  is  no  greater  than  that  of  3,3-dinitrooxetane,  the 
monomer  will  not  undergo  initiation  for  cationic  polymerization.  This  theo¬ 
retical  screening  technique  provides  a  valuable  adjunct  to  the  overall  experi¬ 
mental  program  since  theoretical  calculations  can  be  done  prior  to  synthe¬ 
sis  of  the  monomer.  Many  of  these  monomers  are  very  difficult  to  synthe¬ 
size  and  would  not  be  synthesized  if  there  were  even  reasonable  indication 
that  they  would  not  undergo  initiation  of  cationic  polymerization. 

These  electrostatic  molecular  potential  contour  maps  can  also  be 
used  to  help  answer  question  4)  as  to  good  candidate  partners  for  co¬ 
polymerization.  We  hypothesized  that  since  the  rate  of  polymerization 
is  very  dependent  on  the  basicity  of  the  cyclic  ethers,  for  good  co¬ 
polymerization  the  copolymer  molecules  should  have  reasonably  similar 
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basicities  (as  well  as  ring  strains),  because  if  one  compound  were  much 

more  basic  than  the  other,  it  would  self-polymerize  so  fast  that  it  would 

not  have  much  tendency  to  copolymerize.  Such  copolymers  are  especially 

valuable  in  the  cases  where  the  polymers  formed  from  a  single  monomer 

are  too  brittle  indicating  too  great  a  degree  of  structural  regularity 

18 

(such  as  for  polymers  of  3-fluoro-3-nitrooxetane  ). 

To  shed  insight  on  questions  2)  and  3)  we  are  also  carrying  out  ab- 
initio  MODPOT/VRDDO  and  ab-initio  MODPOT/VRDDO/MERGE  calculations  on  BF3 
approaching  various  series  of  substituted  cyclic  ethers  and  generating 
the  electrostatic  molecular  potential  contour  maps.  These  results  indi¬ 
cate  the  charge  redistribution  and  changes  in  the  maps  as  the  reactants 
approach. 

(a)  3-azidooxetane  (Figures  5  and  6) 

(b)  3-azidooxetane-BF_  complex 

J  o 

B-0  distance  5.0  A  (Figures  6a  and  b) 

(c)  3-azidooxetane-BF,  complex 

B-0  distance  2.  5  A  (Figures  7a  and  b) 

Similar  calculations  are  now  being  carried  out  for  attack  of  H+  on  the 

substituted  cyclic  ethers  since  there  has  been  speculation  by  experimental 

polymer  chemists  that  perhaps  it  is  the  attack  of  H+  (formed  from  a 

Lewis  acid-i^O  complex)  that  triggers  the  polymerization^.  However, 

experimentally  it  has  been  observed  that  different  Lewis  acid  catalysts 

have  different  effects  on  polymerization  of  the  same  monomer.  For  example, 

it  is  possible  to  polymerize  certain  of  these  energetic  monomers  with 

1 9 

PFg  but  not  with  BF^  .  This  would  seem  to  indicate  that  a  catalyst- 
monomer  complex  might  still  be  important.  We  have  suggested  several 
experiment  that  could  help  clarify  the  picture.  First,  to  measure  the 


9 


differing  amounts  of  H+  formed  by  various  Lewis-acid  catalysts  in  the 
polymerization  solvent  mixture  (without  any  monomer  as  well  as  with  mono¬ 
mers)  to  see  if  the  amount  of  H+  formed  could  explain  the  difference  in 
propensity  of  a  monomer  to  polymerize  with  different  Lewis-acid  catalysts. 
Second,  to  add  various  amounts  of  H+  to  solutions  of  monomer  plus  ineffec¬ 
tive  Lewis -acid  catalyst  to  see  if  the  extra  H+  could  trigger  the  reaction. 
Third,  to  add  various  amounts  of  H+  to  solutions  of  monomer  plus 
effective  Lewis-acid  catalysts  to  see  if  the  rate  or  character  of  the 
polymerizations  are  affected. 
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LIST  OF  FIGURE  CAPTIONS 


Figure  1.  Sketches:  Oxetane,  3-Fluoro-3-nitrooxetane,  3,3-Dinitrooxetane 
and  3-Azidooxetane 

Figure  2a.  Oxetane  -  Electrostatic  Molecular  Potential  Contour  Map, 
Parallel  to  the  Mean  Molecular  Plane  xz  (In  Plane) 

b.  Oxetane  -  Electrostatic  Molecular  Potential  Contour  Map, 
Parallel  to  the  xy  Plane  (Through  the  0  Atom) 

Figure  3a.  3-Fluoro-3-nitrooxetane  -  Electrostatic  Molecular  Potential 

Contour  Map,  Parallel  to  the  Mean  Molecular  Plane  xz  (In  Plane) 

b.  3-Fluoro-3-nitrooxetane  -  Electrostatic  Molecular  Potential 
Contour  Map,  Parallel  to  the  xy  Plane  (Through  the  0  Atom) 

Figure  4a.  3,3-Dinitrooxetane  -  Electrostatic  Molecular  Potential 
Contour  Map,  Parallel  to  the  Mean  Molecular  Plane  xz 
(In  Plane) 

b.  3,3-Dinitrooxetane  -  Electrostatic  Molecular  Potential 
Contour  Map,  Parallel  to  the  xy  Plane  (Through  the  0  Atom) 

Figure  5a.  3-Azidooxetane  -  Electrostatic  Molecular  Potential  Contour 
Map,  Parallel  to  the  Mean  Molecular  Plane  xz  (In  Plane) 

b.  3-Azidooxetane  -  Electrostatic  Molecular  Potential  Contour 
Map,  Parallel  to  the  xy  Plane  (Through  the  0  Atom) 

0 

Figure  6a.  3-Azidooxetane-BF3  Complex,  B-0  Distance  5.0  A  -  Electrostatic 
Molecular  Potential  Contour  Map,  Parallel  to  the  Mean  Molecular 
Plane  xz  (In  Plane) 

0 

b.  3-Azidooxetane-BF3  Complex,  B-0  Distance  5.0  A  -  Electrostatic 
Molecular  Potential  Contour  Map,  Parallel  to  the  xy  Plane 
(Through  the  0  Atom) 

o 

Figure  7a.  3-A2idooxetane-BF3  Complex,  B-0  Distance  2.5  A  -  Electrostatic 
Molecular  Potential  Contour  Map,  Parallel  to  the  Mean  Molecular 
Plane  xz  (In  Plane) 

0 

b.  3-Azidooxetane-BF3  Complex,  B-0  Distance  2.5  A  -  Electrostatic 
Molecular  Potential  Contour  Map,  Parallel  to  the  xy  Plane 
(Through  the  0  Atom) 


/ 


I 


Figure  2a  - 


/ 


Figure  6a  - 


Figure  6b  - 


